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^ , Abstract 

C ■ A mean field calculation is carried out to obtain the equation of state (EoS) of 

nuclear matter from a density dependent M3Y interaction (DDM3Y). The energy 
(N ■ per nucleon is minimized to obtain ground state of the symmetric nuclear matter 

^ ■ (SNM). The constants of density dependence of the effective interaction are obtained 

, by reproducing the saturation energy per nucleon and the saturation density of 

\ SNM. The energy variation of the exchange potential is treated properly in the 

negative energy domain of nuclear matter. The EoS of SNM, thus obtained, is not 
only free from the superluminosity problem but also provides excellent estimate of 
\ nuclear incompressibility. The EoS of asymmetric nuclear matter is calculated by 

' adding to the isoscalar part, the isovector component of M3Y interaction. The SNM 

and pure neutron matter EoS are used to calculate the nuclear symmetry energy 
^ , which is found to be consistent with that extracted from the isospin diffusion in 

^ ' heavy-ion collisions at intermediate energies. The /3 equilibrium proton fraction 

calculated from the symmetry energy and related theoretical findings are consistent 
with the constraints derived from the observations on compact stars. 
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1 Introduction 



The investigation of constraints for the high baryonic density behaviour of 
nuclear matter has recently received new impetus with the plans to construct 
a new accelerator facility (FAIR) at GSI Dramstadt. The stiffness of a nu- 
clear equation of state (EoS) is characterised by nuclear incompressibility [1] 
which can be extracted experimentally. Nuclear incompressibility [2,3] also 
determines the velocity of sound in nuclear medium for predictions of shock 
wave generation and propagation. The EoS is of fundamental importance in 
the theories of nucleus-nucleus collisions at energies where the nuclear in- 
compressibility Kq comes into play as well as in the theories of supernova 
explosions [4]. A widely used experimental method is the determination of 
the nuclear incompressibility from the observed giant monopole resonances 
(GMR) [5,6]. Other recent experimental determinations are based upon the 
production of hard photons in heavy ion collisions [7] and from isoscalar giant 
dipole resonances (ISGDR) [8,9,10]. From the experimental data of isoscalar 
giant monopole resonance (ISGMR) conclusion can be drawn that Kq a; 240 
± 20 MeV [11]. The general theoretical observation is that the non-relativistic 
[12] and the relativistic [13] mean field models [14] predict for the bulk incom- 
pressibility for the SNM, Kq, vahies which are significantly different from one 
another, viz. ~ 220-235 MeV and ~ 250-270 MeV respectively. Theoretical 
EoS for the SNM that predict higher Kq values ~ 300 MeV are often called 
'stiff' EoS whereas those EoS which predict smaller Kq values ~ 200 MeV are 
termed as 'soft' EoS. 

The nuclear symmetry energy (NSE) is an important quantity in the equa- 
tion of state of isospin asymmetric nuclear matter. This currently unknown 
quantity plays a key role to the understanding of the structure of systems as 
diverse as the neutron rich nuclei and neutron stars and it enters as an input 
to the heavy ion reactions [15,16]. In general, they can be broadly classified 
into two different forms. One, where the NSE increases monotonically with 
increasing density ('stiff' dependence) [17] and the other, where the NSE in- 
creases initially up to normal nuclear density or beyond and then decreases at 
higher densities ('soft dependence') [18]. Determination of the exact form of 
the density dependence of the NSE is rather important for studying the struc- 
ture of neutron rich nuclei, and studies relevant to astrophysical problems, 
such as the structure of neutron stars and the dynamics of supernova collapse 
[19]. 'Stiff' density dependence of the NSE is predicted to lead a large neu- 
tron skin thickness compared to a 'soft' dependence and can result in rapid 
cooling of a neutron star, and a larger neutron star radius, compared to a 
soft density dependence. A somewhat 'stiff' EoS of nuclear matter need not 
lead to a 'stiff' density dependence of the NSE. Modern constraints from mass 
and mass-radius-relation measurements require a stiff EoS at high densities 
whereas flow data from heavy-ion collisions seem to disfavour too stiff behavior 
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of the EoS. As a compromise hybrid EoS [20,21,22] with a smooth transition 
at high density to quark matter are being proposed. 

In view of rather large differences between the various calculations of the NSE 
present even at the subsaturation densities, the question arises whether one 
can obtain empirical constraints from finite nuclei. Since the degree of isospin 
diffusion in heavy-ion collisions at intermediate energies is affected by the stiff- 
ness of the NSE, these reactions, therefore, can also provide constraints on the 
low energy behaviour of the NSE [17]. However, the high density behaviour 
remains largely undetermined since hardly any data on the simultaneous mea- 
surements of masses and corresponding radii of neutron stars exist whereas 
they can be obtained theoretically by solving Tolman-Oppenheimcr-Volkov 
equation. However there exist indirect indications such as the neutron star 
cooling process. Recently search for the experimental signatures of the moder- 
ately high density behaviour of the NSE has been proposed [18] theoretically 
using several sensitive probes suh as the 7i~ to vr"*" ratio, tranverse collective 
flow and its excitation function as well as the neutron-proton differential flow. 

In the present work, we apply recently discovered astrophysical bounds of high 
density behaviour of nuclear matter from compact star cooling phenomenology 
and also show that the theoretical description of nuclear matter using the den- 
sity dependent M3Y-Reid-Elliott effective interaction [23,24] gives a value of 
nuclear incompressibility which is in excellent agreement with values extracted 
from experiments. The velocity of sound does not become superluminous since 
the energy dependence is treated properly for the negative energy domain of 
nuclear matter. It also provides a symmetry energy that is consistent with 
the empirical value extracted from the measured masses of nuclei. The micro- 
scopic proton-nucleus interaction potential is obtained by folding the density 
of the nucleus with DDM3Y effective interaction whose density dependence 
is determined completely from the nuclear matter calculations. The quantum 
mechanical tunneling probability is calculated within the WKB framework us- 
ing these nuclear potentials. These calculations provide reasonable estimates 
for the observed proton radioactivity lifetimes. Along with earlier works using 
the same formalism, present work provides a unified description of radioactiv- 
ity, scattering, nuclear EoS and NSE. 



2 The nucleeir equation of state for symmetric nuclear matter 

In the present work, density dependence of the effective interaction, DDM3Y, 
is completely determined from nuclear matter calculations. The equilibrium 
density of the nuclear matter is determined by minimizing the energy per 
nucleon. In contrast to our earlier calculations for the nuclear EoS where the 
energy dependence of the zero range potential was treated as fixed at a value 
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corresponding to the equilibrium energy per nucleon eo [25] and assumed to 
vary negligibly with e inside nuclear fluid, in the present calculations the energy 
variation of the zero range potential is treated more accurately by allowing 
it to vary freely but only with the kinetic energy part e'^*" of the energy per 
nucleon e over the entire range of e. This is not only more plausible, but also 
yields excellent result for the incompressibility Kq of the SNM which no more 
suffers from the superluminosity problem. 

The constants of density dependence are determined by reproducing the sat- 
uration conditions. It is worthwhile to mention here that due to attractive 

character of the M3Y forces the saturation condition for cold nuclear matter 
is not fulflllcd. However, the realistic description of nuclear matter proper- 
ties can be obtained with this density dependent M3Y effective interaction. 
Therefore, the constants of density dependence have been obtained by repro- 
ducing the saturation energy per nucleon and the saturation nucleonic density 
of the cold SNM. Based on the Hartrcc or mean field assumption and using the 
DDM3Y interaction, the expression for the energy per nucleon for symmetric 
nuclear matter e is given by 

^ ^ 10m ^ ^ 2 ^ ^ ' 



where Fermi momentum kp — (1.57r^p)5, m is the nucleonic mass equal to 
938.91897 McV/c^ and J^qo represents the volume integral of the isoscalar 
part of the M3Y interaction supplemented by the zero-range potential having 
the form 



= 7999— - 2134-— + Joo(l - ae'^'") where Jqo = -276 MeV fm^ (2) 

where e*^*" = is the kinetic energy part of the energy per nucleon e given 
by Eq.(l). 

The isoscalar t^^^ and the isovector t^^^ components of M3Y interaction po- 
tentials [24,26] supplemented by zero range potentials are given by t^Q^^ {s, e) — 
799922PtzM_2l3422PgM_ 276(1 -ae)5(s) andt^^y^^^^) ^ _4886^^ei_M + 
1176^^2gM + 228(1 - ae)5{s) respectively, where the energy dependence 
parameter «=0.005/MeV. The DDM3Y effective NN interaction is given by 
Voi(s,p, e) = t^^^ {s,e)g{p) where the density dependence g{p) = C(l — /3p") 
and the constants C and /3 of density dependence have been obtained from 
the saturation condition |^ = at p = po and e = eo where po and eo are the 
saturation density and the saturation energy per nucleon respectively. Eq.(l) 
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can be differentiated with respect to p to yield equation 



The equihbrium density of the cold SNM is determined from the saturation 
condition. Then Eq.(l) and Eq.(3) with the saturation condition |^ = can 
be solved simiiltancously for fixed values of the saturation energy per nucleon 
Co and the saturation density po of the cold SNM to obtain the values of j3 and 
C . The constants of density dependence /? and C, thus obtained, are given by 

3 [(1 - + - 7)]A7" 

" [(3n + 1) - (n + l)p + - f )] 



[lOmeo] 2aeoJoo 



where J°qo = ^,;oo(£o*") which means Jyoo at e''"'" = eQ*", the kinetic energy 
part of the saturation energy per nucleon of SNM, kpo — [I.Stt^po]^^^ and 



e^!Mb^ (6) 

5mJOooPo[l -{n + 1)M - ' 



respectively. It is quite obvious that the constants of density dependence C 

and P obtained by this method depend on the saturation energy per nucleon 
Co, the saturation density po, the index n of the density dependent part and 
on the strengths of the M3Y interaction through the volume integral J°oo- 



3 The incompressibility of symmetric nuclear matter 



The incompressibility or the compression modulus of symmetric nuclear mat- 
ter, which is a measure of the curvature of an EoS at saturation density and 
defined as \kF=kF > measures the stiffness of an EoS. The |4 is given by 



dp^ 



Ibmp^ 



-aJooC[l-(n + l)/3p-][-^] + [ 

omp 



^aJooC(l - I3p'")h'kl 
SOmp 



(7) 
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and therefore the incompressibihty Kq of the cold SNM which is defined as 

2 n 2 



can be theoretically obtained as 



9J^,,Cn{n + l)Pp^ 



n+l 



5m 2 
-9aJooC[l-{n + l)Pp^][ 



5m 



+ [ 



3poaJooC{l 



10m 



(9) 



The calculations are performed using the values of the saturation density 
Po=0.1533 fm~^ [2] and the saturation energy per nucleon eo=-15.26 MeV 
[27,28] for the SNM obtained from the co-efficient of the volume term of Bethe- 
Weizsacker mass formula [29,30] which is evaluated by fitting the recent ex- 
perimental and estimated atomic mass excesses from Audi-Wapstra-Thibault 
atomic mass table [31] by minimizing the mean square deviation incorporating 
correction for the electronic binding energy [32]. In a similar recent work, in- 
cluding surface symmetry energy term, Wigner term, shell correction and pro- 
ton form factor correction to Coulomb energy also, turns out to be 15.4496 
MeV [33] (a^ =14.8497 MeV when and A^/^ terms are also included). Using 
the usual values of q;=0.005 MeV~^ for the parameter of energy dependence 
of the zero range potential and n—2/3, the values obtained for the constants 
of density dependence C and f3 and the SNM incompressibihty are 2.2497, 
1.5934 fm^ and 274.7 MeV respectively. The saturation energy per nucleon is 
the volume energy coefficient and the value of -15.26±0.52 MeV covers, more 
or less, the entire range of values obtained for for which now the values 
of C=2.2497±0.0420, /3=1.5934±0.0085 fm^ and the SNM incompressibihty 
Xo=274.7±7.4 MeV. 

The theoretical estimate Kq of the incompressibihty of infinite SNM obtained 
from present approach using DDM3Y is about 270 MeV. The theoretical es- 
timate of Kq from the refractive a-nucleus scattering is about 240-270 MeV 
[34,35] and that by infinite nuclear matter model (INM) [36] claims a well de- 
fined and stable value of Kq = 288 ± 20 MeV and present theoretical estimate 
is in reasonably close agreement with the value obtained by INM which rules 
out any values lower than 200 MeV. Present estimate for the incompressibil- 
ity Kq of the infinite SNM is in good agreement with the experimental value 
of Kq = 300 ± 25 MeV obtained from the giant monopole resonance (GMR) 
[5] and with the the recent experimental determination of Kq based upon the 
production of hard photons in heavy ion collisions which led to the experimen- 
tal estimate of Kq — 290 ± 50 MeV [7] . However, the experimental values of 
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Ko extracted from the isoscalar giant dipolc resonance (ISGDR) arc claimed 
to be smaller [10]. The present non-rclativistic mean field model estimate for 
the nuclear incompressibility Kq for SNM using DDM3Y interaction is rather 
close to the theoretical estimates obtained using relativistic mean field models 
and close to the lower hmit of the older experimental values [5] and close to 
the upper limit of the recent values [6] extracted from experiments. 



Considering the status of experimental determination of the SNM incompress- 
ibihty from data on the compression modes ISGMR and ISGDR of nuclei it 
can be inferred [11] that due to violations of self consistency in HF-RPA 
calculations of the strength functions of giant resonances result in shifts in 
the calculated values of the centroid energies which may be larger in magni- 
tude than the current experimental uncertainties. In fact, the prediction of 
Kq lying in the range of 210-220 MeV were due to the use of a not fully 
self-consistent Skyrme calculations [11]. Correcting for this drawback, Skyrme 
parmetrizations of SLy4 type predict Kq values in the range of 230-240 MeV 
[11]. Moreover, it is possible to build bona fide Skyrme forces so that the 
SNM incompressibility is close to the relativistic value, namely 250-270 MeV. 
Therefore, from the ISGMR experimental data the conclusion can be drawn 
that Ko ^ 240 ± 20 MeV. The ISGDR data tend to point to lower values 
[8,9,10] for Kq. However, there is consensus that the extraction of Kq is in this 
case more problematic for various reasons. In particular, the maximum cross- 
section for ISGDR decreases very strongly at high excitation energy and may 
drop below the current experimental sensitivity for excitation energies [11] 
above 30 and 26 MeV for ^^^Sn and ^°^Pb, respectively. The present value of 
274.7±7.4 MeV for the incompressibility Kq of SNM obtained using DDM3Y 
interaction is, therefore, an excellent theoretical result. 



The constant of density dependence /5=1.5934±0.0085 fm^, which has the di- 
mension of cross section for n—2/3, can be interpreted as the isospin averaged 
effective nucleon-nucleon interaction cross section in ground state symmetric 
nuclear medium. For a nucleon in ground state nuclear matter kp ^ 1.3 fm^^ 
and go ~ hkpc ^ 260 MeV and the present result for the 'in medium' effec- 
tive cross section is reasonably close to the value obtained from a rigorous 
Dirac-Brueckner-Hartree-Fock [37] calculations corresponding to such kp and 
go values which is ~ 12 mb. Using the value of the constant of density de- 
pendence /?=1.5934±0.0085 fm^ corresponding to the standard value of the 
parameter n=2/3 along with the nucleonic density of 0.1533 fm^^, the value 
obtained for the nuclear mean free path A is about 4 fm which is in excellent 
agreement [38] with that obtained using another method. 
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4 The nucleeir equation of state for asymmetric nuclear matter 



The EoS for asymmetric nuclear matter is calculated by adding to the isoscalar 
part, the isovector component [39] of M3Y interaction [40] that do not con- 
tribute to the EoS of SNM. The EoS for SNM and pure neutron matter (PNM) 
are then used to calculate NSE. In this section, implications for the density 
dependence of this NSE in case of neutron stars are discussed, and also pos- 
sible constraints on the density dependence obtained from finite nuclei are 
compared. 

Assuming interacting Fermi gas of neutrons and protons, with isospin asym- 
metry X = y^^, p = Pn + Pp) where p„, pp and p are the neutron, 
proton and nucleonic densities respectively, the energy per nucleon for isospin 
asymmetric nuclear matter can be derived as 



<P.X)^[^\F{X) + {P^){l-Pp-) 



l^Mx) - "i^^. ] (10) 

lOmeo 



where kp — (I.Stt^p) s which is equal to Fermi momentum in case of SNM, the 
kinetic energy per nucleon e*^*" = [^^]F{X) with F{X) = [(l±^)!^!±ilz^lf^!] 
and Jy = Jvoo -\- Jvoi and J^oi represents the volume integral of the isovector 
part of the M3Y interaction supplemented by the zero-range potential having 
the form 



ATT 4:77 

= -4886— + 1176— r + Joi(l - ae'^'") where Jqi = 228 MeV fm^ (11) 
4'' 2.5'^ 

where J — Joo + X^Joi and 
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The pressure P and the energy density s of nuclear matter with isospin asym- 
metry X are given by 



P 



. de 



[p^]F{X)+p^-^[l-{n + lW^] 



dp '■'^ 5m 
-p^a JC[1 - 



Win 



(14) 



and 



e = p{e + mc^) 



Pii^-^MX) + (^)(l - Ppn + mc% (15) 



respectively, and thus the velocity of sound Vg in nuclear matter with isospin 
asymmetry X is given by 



c 




[2P| 



[e + mc2 + p|] 



(16) 



The incomprcssibilities for isospin asymmetric nuclear matter arc evaluated 
at saturation densities ps with the condition |^ = which corresponds to 
vanishing pressure. The incompressibility Kg for isospin asymmetric nuclear 
matter is therefore expressed as 



' ^ 5^ ^^ 2 ^ ' 

om lUm 

Here kp^ means that the kp is evaluated at saturation density ps- = J'Iqq + 
J*'q^ is the J,^ at e'^*" = ej*" which is the kinetic energy part of the saturation 
energy per nucleon and J = Jqo + X'^J^ 



01- 



In Table- 1 incompressibility of isospin asymmetric nuclear matter Kg as a 
function of the isospin asymmetry parameter X, using the usual value of 
n=2/3 and energy dependence parameter a=0.005 MeV^^, is provided. The 
magnitude of the incompressibility Kg decreases with the isospin asymmetry 
X due to lowering of the saturation densities pg with X as well as decrease 
in the EoS curvature. At high isospin asymmetry X, the isospin asymmetric 
nuclear matter does not have a minimum signifying that it can never be bound 
by itself due to nuclear interaction. However, the j3 equlibrated nuclear matter 
which is a highly neutron rich asymmetric nuclear matter exists in the core of 
the neutron stars since its E/A is lower than that of SNM at high densities and 
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Table 1 

Incompressibility of isospin asymmetric nuclear matter using the usual value of 
n=2/3 and energy dependence parameter a=0.005 MeV~^. 



X 


Ps 


Ks 




fm~^ 


MeV 


0.0 


0.1533 


274.7 


0.1 


0.1525 


270.4 


0.2 


0.1500 


257.7 


0.3 


0.1457 


236.6 


0.4 


0.1392 


207.6 


0.5 


0.1300 


171.2 



is unbound by the nuclear force but can be bound due to high gravitational 
field realizable inside neutron stars. 

The theoretical estimates of the pressure P and velocity of sound Vg of SNM 
and isospin asymmetric nuclear matter including PNM are calculated as func- 
tions of nucleonic density p and energy density e using the usual value of 
0.005 MeV~^ for the parameter a of energy dependence of the zero range 
potential and also the standard value of the parameter n=2/3. Unlike other 
non-relativistic EoS, present EoS does not suffer from superluminosity at all 
for SNM and for PNM problem of super luminosity occurs only at very high 
densities (p > 12po), higher than those encountered at the centres of neutron 
stars. In our earlier calculations for the nuclear EoS where the energy depen- 
dence of the zero range potential was treated as fixed at a value corresponding 
to the equilibrium energy per nucleon eo [28] and assumed to vary negligibly 
with e inside nuclear fluid caused superluminosity problems for both SNM 
and PNM and at much lower densities [25] hke the EoS obtained using the 
+ TNI interaction [41]. In the present calculations the energy variation of 
the zero range potential is treated more accurately allowing it to vary freely 
but only with the kinetic energy part e^*" of the energy per nucleon e over the 
entire range of e. 

In Fig.-l the energy per nucleon e of SNM and PNM arc plotted as functions of 
p. The continuous lines represent the curves for the present calculations using 
saturation energy per nucleon of -15.26 MeV whereas the dotted lines represent 
the same using V14 + TNI interaction [41] and the dash-dotted hues represent 
the same for the A18 model using variational chain summation (VCS) [42] for 
the SNM and PNM. The minimum of the energy per nucleon equaling the 
saturation energy -15.26 MeV for the present calculations occurs precisely at 
the saturation density po=0.1533 fm~^ since equilibrium density po of the cold 
SNM is determined from the saturation condition f^=0 at p—po and e—eQ. 
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Fig. -2 presents the plots of the energy per nucleon e of nuclear matter with 
different isospin asymmetry X as functions of p/po ioi the present calculations. 
The pressure P of SNM and PNM are plotted in Fig.-3 and Fig.-4 respectively 
as functions of p/po- The continuous lines represent the present calculations 
whereas the dotted lines represent the same using the A18 model using vari- 
ational chain summation (VCS) of Akmal et al. [42] for the SNM and PNM. 
The dash-dotted line of Fig. -3 represents plot of P versus p/po for SNM for 
RMF using NL3 parameter set [43] and the area enclosed by the coninuous 
line corresponds to the region of pressures consistent with the experimental 
flow data [44]. It is interesting to note that the RMF-NL3 incompressibility for 
SNM is 271.76 MeV [45] which about the same as 274.7±7.4 MeV obtained for 
the present calculation. The areas enclosed by the continuous and the dashed 
lines in Fig.-4 correspond to the pressure regions for neutron matter consistent 
with the experimental flow data after inclusion of the pressures from asym- 
metry terms with weak (soft NM) and strong (stiff NM) density dependences, 
respectively [44]. In Fig.-5 the velocity of sound Vg in SNM and PNM and the 
energy density e of SNM and PNM for the present calculations are plotted as 
functions of p/ po. The continuous lines represent the velocity of sound in units 
of 10~^c whereas the dotted lines represent energy density in MeV fm~^. 



5 The nuclecir symmetry energy 



The nuclear symmetry energy Esym{p) represents a penalty levied on the sys- 
tem as it departs from the symmetric limit of equal number of protons and 
neutrons and can be defined as the energy required per nucleon to change the 
SNM to pure neutron matter (PNM) [46] 

^.2/m(p)=e(p,l)-e(p,0) (18) 

and therefore using Eq.(lO) for X = 1 and Eq.(l), the NSE can be given by 
EsyM = (2^/^ - l)^EU^f' + §P(1 - Pp^Jvoi (19) 

O Po z 

where the Fermi energy E'^ — for the SNM at ground state. 
5.1 Nuclear symmetry energy at high baryonic density 

The first term of the right hand side is the kinetic energy contribution with 
density dependence of Oop^^^ whereas the second term arising due to nuclear 
interaction has a density dependence of the form of Oip -|- a2p""^^ + a^p^/^ -\- 
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a^p'^'^^^^ since J,,,oi is a function of e'^*" which varies as p^/^ and 00,01,02,03 
and 04 are constants with respect to the nucleonic density p or the parameter 
n. If one uses an ahernative definition [47] of Egymip) = \^~^x^ \x={) for the 
nuclear symmetry energy, only the term (2^/^ — 1) of the above equation gets 
replaced by 5/9 [which is about five percent less than (2^/^ — 1)] and reduces 
only the kinetic energy contribution. 

In Fig.-6 plots of E/A for SNM, PNM and NSE as functions of p/po are 
shown for n=2/3. The density dependence of the NSE at subnormal density 
from isospin diffusion [17] in heavy- ion collisions at intermediate energies has 
an approximate form of 31. 6[^]^'''^ MeV. This low energy behaviour of NSE 
K, 31.6[^]^ °^ MeV is close to that obtained using Eq.(19) at subnormal den- 
sities. At higher densities the present NSE using DDM3Y interaction peaks 
at p ~ 1.95po and becomes negative at p ~ 4.7po- A negative NSE at high 
densities implies that the pure neutron matter becomes the most stable state. 
Consequently, pure neutron matter exists near the core of the neutron stars. 



5.2 Nuclear symmetry energy at low baryonic density 



The volume symmetry energy coefficient Sy extracted from the masses of fi- 
nite nuclei provides a constraint on the nuclear symmetry energy at nuclear 
density Esym{po). The value of 5"^ = 30.048 ± 0.004 MeV recently extracted 
[48] from the measured atomic mass excesses of 2228 nuclei is reasonably 
close to the theoretical estimate of the value of NSE at the saturation den- 
sity £^syr„(po)=30.71±0.26 MeV obtained from the present calculations us- 
ing DDM3Y interaction. Instead of Eq.(18), if an alternative definition [47] 
lt'2 \x=o of the nuclear symmetry energy is used, then its 
value is 30.03±0.26 MeV. In ref. [49] it is between 29.10 MeV to 32.67 MeV 
and that obtained by the liquid droplet model calculation of ref. [50] is 27.3 
MeV whereas in ref. [51] it is 28.0 MeV. It should be mentioned that the value 
of the volume symmetry parameter Sy in some advanced mass description [52] 
is close to the present value which with their — vol vol ~ Sv equals 29.3 MeV. 
The value of NSE at nuclear saturation density ^ 30 MeV, therefore, seems 
well established empirically. Theoretically different parametrizations of the 
relativistic mean-field (RMF) models, which fit obseravables for isospin sym- 
metric nuclei well, lead to a relatively wide range of predictions 24-40 MeV for 
EgymiPo)- The present result of 30.71±0.26 MeV of the mean field calculation 
is close to the results of the calculation using Skyrme interaction SkMP (29.9 
MeV) [53], Avl8+5v+mX* variational calculation (30.1 MeV) [42] and field 
theoretical calculation DD-F (31.6 MeV) [46]. 
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5.3 Compact star constraints 



The knowledge of the density dependence of nuclear symmetry energy is im- 
portant for understanding not only the stucture of radioactive nuclei but also 
many important issues in nuclear astrophysics, such as nucleosynthesis during 
presupernova evolution of massive stars and the cooling of protoneutron stars. 
A neutron star without neutrino trappings can be considered as a.n,p,e mat- 
ter consisting of neutrons (n), protons (p) and electrons (e). The neutrinos do 
not accumulate in neutron stars because of its very small interaction probabil- 
ity and correspondingly very high mean free path [54,55]. The /3 equilibrium 
proton fraction xp [— Pp/ p] is determined by [56] 

hc{'iT:''pxpf/^ = AEsym{p){l - 2x^). (20) 



The equilibrium proton fraction is therefore entirely determined by the NSE. 
The 13 equilibrium proton fraction calculated using the present NSE is plotted 
as function of p/po in Fig.-7. The maximum of ^ 0.044 occurs at p 
1.35po and goes to zero at p ^ 4.5po for n—2/3. The NSE extracted from the 
isospin diffusion in the intermediate energy heavy-ion collisions, having the 
approximate form of 31. 6[^]^ '^^ McV. provides a monotonically increasing /3 
equilibrium proton fraction and therefore can not be extended beyond normal 
nuclear matter densities. Present calculation, using NSE given by Eq.(19), of 
the P equihbrium proton fraction forbids the direct URCA process since the 
equilibrium proton fraction is always less than 1/9 [56] which is consistent 
with the fact that there are no strong indications that fast cooling occurs. 
Moreover, recently it has been concluded theoretically that an acceptable EoS 
of asymmetric nuclear matter (such as P equilibrated neutron matter) shall not 
allow the direct URCA process to occur in neutron stars with masses below 1.5 
solar masses [46]. Although a recent experimental observation suggests high 
heat conductivity and enhanced core cooling process indicating the enhanced 
level of neutrino emission but that can be via the direct URCA process or 
Cooper-pairing [57]. Also observations of massive compact stars in the mass 
range of 2.1±0.2 solar mass to a la confidence level (and 2.llo;5 solar mass 
to a 2(7 confidence level) and 2.0±0.1 solar mass and the lower bound for the 
mass-radius relation of isolated pulsar RX J1856 imply a rather 'stiff' nuclear 
EoS [46]. The present NSE is 'soft' because it increases initially with nucleonic 
density up to about two times the normal nuclear density and then decreases 
monotonically at higher densities. It is interesting to observe that although 
the SNM incompressibility is slightly on the higher side and the present EoS 
is 'stiff', yet the present calculations provide a rather 'soft' nuclear symmetry 
energy and thus satisfy the astrophysical constraints. 
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Table 2 

Comparison between experimentally measured and theoretically calculated half- 
lives of spherical proton emitters. The asterisk symbol (*) in the experimental Q 
values denotes the isomeric state. The experimental Q values, half lives and I values 
are from ref. [58]. The results of the present calculations using the isoscalar and 
isovector components of DDM3Y folded potentials are compared with the experi- 
mental values and with the results of UFM estimates [59] . Experimental errors in Q 
[58] values and corresponding errors in calculated half-lives are inside parentheses. 



Parent 


I 


Q 


1^* tpt 


2"^^ tpt 


3'"'^ tpt 


Expt. 


This work 


UFM 




n 


MeV 


i?i[fm] 


Ra[im] 


R,[im] 


logioT{s) 


logioT{s) 


logioT{s) 




2 


0.491(15) 


1.43 


6.69 


134.30 


2.049;^^? 


1.90(45) 


2.085 




5 


1.753(10) 


3.20 


6.63 


56.27 




-5.28(7) 


-5.170 




5 


1.071(3) 


3.18 


6.63 


88.65 


o.59il°:li 


0.83(4) 


1.095 




2 


1.139(5) 


1.44 


7.28 


78.97 


-S 444+0-046 


-3.46(6) 


-3.199 




5 


1.283(4) 


3.21 


6.67 


78.23 


-i.i»U_o.o64 


-0.74(4) 


-0.859 


150^^* 


2 


1.317(15) 


1.45 


7.33 


71.79 


A r9q+0.620 


-4.46(15) 


-4.556 




5 


1.255(3) 


3.21 


6.69 


78.41 


-0.89610:011 


-0.82(4) 


-0.573 




2 


1.332(10) 


1.46 


7.35 


69.63 


-4.7961°:°^? 


-4.96(10) 


-4.715 




5 


1.791(10) 


3.21 


6.78 


57.83 


-4 021+0-125 


-4.80(7) 


-4.637 




2 


1.028(5) 


1.47 


7.37 


94.18 


-0 6 20+0-082 


-0.47(8) 


-0.461 




5 


1.130(8) 


3.21 


6.76 


90.30 


n Q4Q+0-100 

u.y^y_o 


1.50(10) 


1.446 







0.947(7) 


0.00 


7.55 


98.95 


n coq+O.lSS 
-U.OZO_o.i98 


-0.51(12) 


-0.126 


i60i?e 


2 


1.284(6) 


1.45 


7.43 


77.67 


3 046+0-075 


-3.08(7) 


-3.109 


i6ii?e 





1.214(6) 


0.00 


7.62 


79.33 


q 4q9+0.045 


-3.53(7) 


-3.231 




5 


1.338(7) 


3.22 


6.84 


77.47 


-0.48810:0| 


-0.75(8) 


-0.458 


164j^ 


5 


1.844(9) 


3.20 


6.91 


59.97 


-3.95910:190 


-4.08(6) 


-4.193 


165 j^* 


5 


1.733(7) 


3.21 


6.93 


62.35 


-3.46910:?^^ 


-3.67(5) 


-3.428 


166 jj. 


2 


1.168(8) 


1.47 


7.49 


87.51 


-0.82410:166 


-1.19(10) 


-1.160 


166 J J,* 


5 


1.340(8) 


3.22 


6.91 


80.67 


-0.076l0:i2| 


0.06(9) 


0.021 


167 j-^ 





1.086(6) 


0.00 


7.68 


91.08 


-0.959l0:0i 


-1.35(8) 


-0.943 


167 j^* 


5 


1.261(7) 


3.22 


6.92 


83.82 


0.87510:?^ 


0.54(8) 


0.890 


171 An 





1.469(17) 


0.00 


7.74 


69.09 


-4.77010:185 


-5.10(16) 


-4.794 


171 An* 


5 


1.718(6) 


3.21 


7.01 


64.25 


-2 654+0-054 


-3.19(5) 


-2.917 


i77r/ 





1.180(20) 


0.00 


7.76 


88.25 


+0.191 
-"-•-■^ '^-0.349 


-1.44(26) 


-0.993 




5 


1.986(10) 


3.22 


7.10 


57.43 


o 7+0.095 
a.a'±(_o.l22 


-4.64(6) 


-4.379 







1.624(16) 


0.00 


7.91 


65.71 


A 22Q+0.068 

-4.zzy_o 081 


-5.53(14) 


-5.184 



6 Folding model analyses using effective interaction whose density 
dependence determined from nuclear matter calculation 



Microscopic proton- nucleus interaction potentials are obtained by single fold- 
ing the density of the nucleus with M3Y effective interaction supplemented 
by a zero-range pseudo-potential for exchange along with the density depen- 
dence. Parameters of the density dependence, C=2.2497 and /9=1.5934 fm^, 
obtained here from the nuclear matter calculations assuming kinetic energy 
dependence of zero range potential, are used. 

The half lives of the decays of spherical nuclei away from proton drip line by 
proton emissions are estimated theoretically. The half life of a parent nucleus 
decaying via proton emission is calculated using the WKB barrier penetration 
probability. The WKB method is found quite satisfactory and even better than 
the S-matrix method for calculating half widths of the a decay of superheavy 
elements [60]. For the present calculations, the zero point vibration energies 
used here are given by eqn.(5) of ref. [61] extended to protons and the ex- 
perimental Q values [58] are used. Spherical charge distributions are used for 
Coulomb interaction potentials. The same set of data of rcf. [59] has been used 
for the present calculations using C=2.2497 and /3=1.5934 fm^ and presented 
in Table-2. The agreement of the present calculations with a wide range of 
experimental data for the proton radioactivity lifetimes are reasonable. 

Since the density dependence of the effective projectile-nucleon interaction was 
found to be fairly independent of the projectile [62], as long as the projectile- 
nucleus interaction was amenable to a single-folding prescription, the density 
dependent effects on the nucleon-nuclcon interaction were factorized into a tar- 
get term times a projectile term and used successfully in case of a radioctivity 
of nuclei [63] including superheavies [64,65,66] and the cluster radioactivity 
[63,67]. The calculations were performed for elastic and inelastic scattering 
of protons from nuclei ^^Ne, ^^O, ^°0, ^^O using C=2.07 and /5=1.624 fm^ 
[68,69]. It is needless to say that the present value of /5=1.5934±0.0085 fm^, 
obtained by treating the energy variation of the zero range potential prop- 
erly by allowing it to vary freely with the kinetic energy, which changes by 
about one percent neither changes the shape of the potential significantly nor 
the quality of fit or the values extracted for the nuclear deformations. How- 
ever, since the value of C, which acts as the overall normalisation constant 
for the nuclear potentials, changes by about six percent, causes changes but 
only to the renormalizations required for the potentials. Therefore, it provides 
reasonable description for elastic and inelastic scattering of protons and the 
deformation parameters extracted from these analyses are in good agreement 
with the quadrupole deformations obtained from the available experimental 
B{E1) values [70]. 
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7 Summciry and conclusion 



A mean field calculation is carried out to obtain the equation of state of 
nuclear matter from a density dependent M3Y interaction (DDM3Y). The 
microscopic nuclear potentials are obtained by folding the DDM3Y effective 
interaction with the densities of interacting nuclei. The energy per nucleon is 
minimized to obtain ground state of the symmetric nuclear matter (SNM). 
The constants of density dependence of the effective interaction are obtained 
by reproducing the saturation energy per nucleon and the saturation density 
of SNM. The EoS of asymmetric nuclear matter is calculated by adding to the 
isoscalar part, the isovector component of M3Y interaction. The SNM and 
pure neutron matter EoS are used to calculate the nuclear symmetry energy 
which is found to be consistent with that extracted from the isospin diffusion in 
heavy-ion collisions at intermediate energies. The microscopic proton-nucleus 
interaction potential is obtained by folding the density of the nucleus with 
DDM3Y effective interaction whose density dependence is determined com- 
pletely from the nuclear matter calculations. 



In this work the energy variation of the exchange potential is treated properly 
in the negative energy domain of nuclear matter. The EoS of SNM, thus ob- 
tained, is free from the super luminosity problem encountered in some previous 
prescriptions. Moreover, the result of the present calculation for the compres- 
sion modulus for the infinite symmetric nucler matter is in better agreement 
with that extracted from experiments. The calculated (3 equilibrium proton 
fraction forbids direct URCA process which is consistent with the fact that 
there arc no strong indications that fast cooling occurs. The results of the 
present calculations using single folded microscopic potentials for the proton- 
radioactivity lifetimes are in good agreement over a wide range of experimental 
data. We find that it also provides reasonable description for the elastic and 
inelastic scattering of protons and the deformation parameters extracted from 
the analyses are in good agreement with the available results. The results of 
the present calculations using microscopic potentials for half life calculations of 
a decays are found to be in excellent agreement with experimental data. These 
calculations also provide reliable estimates for the observed a decay lifetimes 
of the newly synthesized superheavy elements. It is, therefore, pertinent to 
conclude that a unified description of the symmetric and asymmetric nuclear 
matter, elastic and inelastic scattering, and cluster, a and proton radioactiv- 
ities is achieved. With the energies and interaction rates foreseen at FAIR, 
the compressed baryonic matter (CBM) will create highest baryon densities 
in nucleus-nucleus collisions to explore the properties of superdense baryonic 
matter and the in-medium modifications of hadrons. 
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400 




p (fm-3) 

Fig. 1. The energy per nucleon e = E/A of SNM (spin and isospin symmetric nuclear 
matter) and PNM (pure neutron matter) as functions of p. The continuous hues 
represent curves for the present calculations using DDM3Y interaction, the dotted 
lines represent the same using vn + TNI interaction [41] and the dash-dotted lines 
represent the same for A18 model using variational chain summation (VCS) [42]. 
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Fig. 2. The energy per nucleon e=E/A of nuclear matter with different isospin asym- 
metry X as functions of p/po for the present calculations using DDM3Y interaction. 
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Fig. 3. The pressure P of SNM (spin and isospin symmetric nuclear matter) as 
a function of p/po- The continuous hnes represent the present calculations using 
eo = -15.26ib0.52 MeV. The dotted line represents the same using the A18 model 
using variational chain summation (VCS) of Akmal et al. [42], the dash-dotted line 
represents the RMF calculations using NL3 parameter set [43] whereas the area 
enclosed by the coninuous line corresponds to the region of pressures consistent 
with the experimental flow data [44] for SNM. 
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1 2 3 4 5 6 7 

P/Po 

Fig. 4. The pressure P of PNM (pure neutron matter) as a function of p/ Pq. The 
continuous hnes represent the present calculations using eo = -15.26±0.52 MeV. 
The dotted line represents the same using the A18 model using variational chain 
summation (VCS) of Akmal et al. [42] whereas the areas enclosed by the continuous 
and the dashed lines correspond to the pressure regions for neutron matter consistent 
with the experimental flow data after inclusion of the pressures from asymmetry 
terms with weak (soft NM) and strong (stiff NM) density dependences, respectively 
[44]. 
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P/Po 

Fig. 5. The velocity of sound Vg in SNM (spin and isospin symmetric nuclear matter) 
and PNM (pure neutron matter) and the energy density e of SNM and PNM as 
functions of p/po for the present calculations using DDM3Y interaction. The con- 
tinuous lines represent the velocity of sound in units of lO^^c whereas the dotted 
lines represent energy density in MeV fm^^. 
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Fig. 6. The energy per nucleon e=E/A of SNM (spin and isospin symmetric nuclear 
matter), PNM (pure neutron matter) and NSE (nuclear symmetry energy Esym) are 
plotted as functions of p/po for the present calculations using DDM3Y interaction. 
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Fig. 7. The /? equilibrium proton fraction calculated with NSE (nuclear symmetry 
energy) obtained using DDM3Y interaction is plotted as a function of p/po- 
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